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Abstract 

We propose that some puzzles in semileptonic decays of B mesons to the broad D 
states could be clarified by studying at LHCb the corresponding decays with strange mesons 
B° — > D~ . In particular, we point out that the non-leptonic decay B° — > D~ tt + and the 
like, being Class I decays (where factorization is expected to hold), could be a first step in 
this direction. The interpretation of results in both semileptonic and non-leptonic decays will 
presumably be easier due to the narrowness of the D~ state. On the other hand, we make 
a careful and detailed study of the experimental and theoretical situation in the case of the 
wide non-strange D case, and we update previous analyses. 



1 Motivation 

The long-standing problem of weak transitions between B and the broad L = 1 (j = 1/2) states 
D** remains interesting to elucidate for at least two reasons: 

1. A lot of theoretical effort has been devoted to understand these transitions by using several 
different approaches. 

1 Note that we use the spectroscopic labels related to the heavy quark limit in which the angular momentum of 
the light degrees of freedom j is a good quantum number. In that limit L — 1 corresponds to both j p = (l/2) + and 
j p = (3/2) + doublets of heavy-light mesons. The (l/2) + -doublet is denoted as [Dq,D[], while the (3/2) + -doublet 
is referred to as [DijD^]- The states with the strange valence quark are distinguished by an extra index "a". 
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2. A considerable experimental effort to measure the corresponding quantities lead to contro- 
versies: the experiments seemed to disagree among themselves and/or with theory. 

In the following, we explain the current situation in Part HI by discussing the theoretical expec- 
tations, experimental results and by comparing theory with experiment. We distinguish between 
the semileptonic and non-leptonic decays. In Part [TT] we propose a way to clarify the puzzles by 
studying the strange -D**-states which happen to be narrow. 

Part I 

Difficulties with B — > D weak decays 



Many papers and notes have been devoted to the above problem. The issues have been discussed 
and summarized some years ago in refs. [U [2] to which we refer for complementary information 
and references. 

2 The broad L = 1 (J = 1/2) c(u,d) D** states 

There are two states with L = 1 (j = 1/2): one with J p = + (Do), the other with J p = 1 + 
(D[). Both are expected to be broad, because of the strong S*-wave decays to D^tt, and the fact 
that their mass is expected to be notably above the D^tt thresholds. 

They have been most clearly observed in the non-leptonic B — >■ D tt decays wherefrom their 
properties, like widths and masses, have been established. Although the semileptonic decay 
rates are much larger than the non-leptonic ones, the number of observed events is in a reversed 
proportion as we explain below. 

The decay rates of the two D** states into D^n are identified as their total widths, which 
is roughly expected from simple quark model calculations [3]. The identification of the very 
broad bumps in D^tt with the expected D** states is plausible although (i) the identification of 
very broad resonances is not safe, and (ii) the observed discrepancies between the predicted and 
observed D** states made the eg- interpretation questionable (q being either u, d or s quark). We 
will briefly return to the latter point in Sec. [5j 

Similar to the above-mentioned bumps were observed in the semileptonic B — >■ D ir£ + vg 
decays, but not with an accuracy allowing to determine the resonance's features independently. 
Rather, one uses the D** properties found in B — > D tt as input in order to estimate the 
semileptonic decay rate. 

Theoretically, however, the semileptonic decays are simpler to describe and require less as- 
sumptions than the non-leptonic ones, and we will discuss them in that order. 

3 Theoretical predictions for B — > D IW functions 

In the heavy quark limit for the c and b quarks all the form factors governing B — > D £ + V£ decays 
are related by simple relations and proportional to one of the two Isgur-Wise (IW) functions, 
T i/2( w ) f° r the final hadron belonging to the j = 1/2 doublet, or t^/2(w) for D** being one of 
the mesons from the j = 3/2 doublet. These functions parameterize the non-perturbative QCD 
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dynamics of the vector or axial current matrix elements [5] as, for example, 
(0+ \A^\ 0") = --r^r(<V " t 1/2 (w) , 

(2 + \A,\ 0-) = ^ -J= [(1 + - ^"^ey r 3/2 H , (1) 



where v, v' are the velocity vectors of the initial and final mesons, 6^ is the polarisation tensor 
of the 2 + -state, and w = v ■ v' . The normalisation of states is (27r) 3 5(v — if). 

The argument of Tj(w) varies between 1 < w J$ 1.3, as it can be easily seen from 



2 tub mD** 

for g^ in = m 2 ps 0, and q^ax = i m B — mn**) 2 . For the non-leptonic decays q 2 = m 2 is fixed and 
corresponds to w « 1.3. Importantly, Ti/ 2 (w) is known to be a slowly varying function of it!, and 
it is a common practice to focus on its normalization at zero recoil w = 1, namely ^^(l). For 
example, in ref. [1] it was found that Ti/ 2 (w) = Ti/ 2 (1)[1 — 0.83(w; — 1) + . . . ]. 

3.1 Inclusive sum rules 

A useful constraint concerning the values of Tj(l) is provided by what we can call Bjorken-like 
or inclusive sum rules, which are not to be confused with the "QCD Sum Rules" a la SVZ [6] 
in that they do not pretend to go beyond equalling the sum over all states of suitable quantum 
numbers to the result obtained by employing the operator product expansion (OPE). They in 
fact reflect the duality with free quarks. One of the most famous such sum rule is the so called 
Uraltsev SR [7], 

Ei-3 ( /1d)r-i^a)i 2 =i> (3) 

with n labeling possible radial excitations (n = being the ground state). Focusing only onto the 
ground states suggests the inequality |ti /2 (1) | < 1 73/2(1) I j which is also confirmed by the similar 
sum rule studied in ref. [S]. This is obviously not a theorem, but relies on assumption that the 
lowest state dominates in each channel. The right hand side, 1/4 may seem a small difference, 
but since |t"i/2(1)| 2 i s a small number, the ratio |T3/2(l)/ r i/ 2 (^)l ^ s ra ^ ner large, 

tagW!! = i + v 4 ( 4) 

|r 1/2 (D| 2 + ln/ 2 (DI 2 ' U 

when considering the lowest states only. This tendency is observed in actual theoretical calcula- 
tions, except in the QCD sum rule calculation of ref. [9]. In the semileptonic rates it is further 
exacerbated by the kinematic factors. 

3.2 Lattice QCD predictions 

The only method allowing to compute these form factors, strictly based on QCD, is the method 
of numerical simulations of QCD on the lattice. The first calculation of 71/2(1) has been made in 
ref. [10] and then extended and improved in ref. [llj , where the computation is made by including 
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the Nf = 2 flavors of dynamical ("sea") quarks. The results of ref. obtained at a single lattice 
spacing, exhibit a negligible dependence on the light quark mass and read: 

r 3/2 (l) = 0.528(23) , r 1/2 (l) = 0.297(26) , (5) 

where the errors do not include the discretization nor the finite volume effects. Note also that 
one cannot easily calculate these form factors away from w = 1 on the lattice. 

3.3 Quark model predictions 

Familiar opinions that "any model would do" or that "you may get anything you want by choos- 
ing a suitable modeT come from disregarding the necessary careful discussions which allow to 
estimate the overall merits of respective models by consideration of the largest possible set of 
phenomenological data and of theoretical consistency and inputs. 

There is no perfect model, other than QCD, but there are definitely bad models and more 
satisfactory ones. One necessary general feature is that for heavy-light systems they should be 
relativistic. As to external motion of hadrons one can use the Bakamjian-Thomas (BT) approach 
which provides a definite way to define states in motion starting from states at rest by constructing 
an explicit Poincare algebra. A particular case is obtained by performing boosts to the infinite 
momentum frame, which gives the familiar null-plane formalism. Covariance of current matrix 
elements is ensured in the heavy mass limit only. Note that the above inclusive sum rules, required 
by QCD, are exactly satisfied by the BT quark model approach. 

Within the BT quark model approach the difference between t 3/ / 2 (1) and T]y 2 (l) comes from 
the Wigner rotations of the light spectator quark, which acts differently for j = 1/2 and for 
j = 3/2 states. One finds that the difference |t 3/ / 2 (1)| — |ti/ 2 (1)| is positive and large [3J. 

In addition to the quark model framework, one also has to choose a (necessarily relativistic) 
potential model to fix the wave functions at rest.H The guiding principle in choosing the potential 
is obviously the requirement to describe as broad range of observed hadrons as possible. In 
that respect, the standard Godfrey-Isgur (GI) potential model provides the best description of 
the whole spectroscopy. By using the wave functions fixed by the GI potential model, the BT 
approach leads to the following results: 

r 3/2 (l)~0.54, r 1/2 (l)~0.22. (6) 

The agreement with the results of lattice calculations ([5]), which have been produced much later, 
is striking. The suppression of Ti/ 2 (1) with respect to t 3 / 2 (1) could be even stronger if other 
potentials (other than GI) are chosen, while t 3 / 2 (1) remains stable. We do not quote errors to 
the above results because there is no clearly admitted definition of errors in the quark models, 
unlike in the well defined method of lattice QCD. For instance, it would not make much sense to 
make an arbitrary variation of parameters without taking into account the whole set of possible 
phenomenological applications, most of which depend on additional modeling. 

Before continuing we would like to emphasize the consistency of the results obtained in the 
static limit of QCD on the lattice with the results obtained by using the BT framework with a 
suitable potential model. Such an agreement is not just a matter of luck. A similar agreement 
has been observed in a very detailed manner for the distribution of the axial, scalar and vector 

2 In a very extensive work, H. Cheng et al. [T5] have made predictions for the transitions to the D**-states in the 
null-plane formalism, including the finite m&, c effects, which is quite useful. However, to be conclusive, a necessary 
step in this approach, which remains to be done, would be to systematically deduce the wave functions from a 
relativistic potential model constrained by the spectrum. 
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charges in the static-light mesons with either L = 0, or L = 1 |13j . The advantage of quark 
models is that one can easily calculate the w-dependence of Ti/ 2j 3/ 2 (u>), needed when computing 
the branching ratios, and get moreover an intuitive insight. 

3.4 QCD sum rules approach to form factors 

The results from QCD sum rules are less safe and less intuitive, and the results for Tj(l) presented 
so far in the literature do not agree among themselves. A major concern is that the results depend 
quite strongly on the choice of the interpolating field for the L>**-states. 
Results of the first calculations presented in refs. (HUH], El 

r 3/2 (l)~0.25 > r 1/2 (l)~ 0.35(8), (7) 

clearly challenge the hierarchy |7i/ 2 (l)| < 1 73/2(1) I- A little later another QCD sum rules compu- 
tation resulted in [15], 

r 3/2 (l) ~ 0.43(8) , r 1/2 (l) ~ 0.13(4) , (8) 

arguing that the usual local scalar interpolating field operator does not lead to a satisfactory 
sum rule, due to a lack of perturbative contribution. To circumvent the problem they used 
the operators with covariant derivative instead. It must be stressed that the quoted "errors" 
in eqs. (|7l8p are not errors in the usual sense of indicating a possible deviation from the true 
value. They merely indicate the variation of the result within the chosen range for the continuum 
threshold. Therefore, one should neither consider 7"i/ 2 (l) — 0.13(4) as being incompatible with 
the result in eq. 0, nor incompatible with the values given in eqs. (|5l6p . The difference between 
the values in eqs. ([7|) and (JSJ) could be viewed as an indicator of a possible uncertainty of the 
method. What is to be actually retained from the results of ref. [15] is that the hierarchy is 
similar to the one found in the lattice QCD and in the quark model discussed above. 

3.5 Phenomenology with 7i/ 2 (l) and r 3 / 2 (l) 

From the above discussion we see that there is a growing evidence that the Uraltsev sum rule is 
well respected by the actual values for the IW functions involving the n = Z)**-states at w = 1, 
and that t 1 / 2 (1) < t 3 / 2 (1). Of course the discussion so far has been restrained to the heavy quark 
limit of QCD. The impact of the corrections arising from the finiteness of the heavy quark mass 
has not been much discussed in the literature and there is no available lattice QCD result that 
would help us assess the size of these corrections. An early careful estimate of these corrections 
within a systematic HQET expansion of ref. [16j suggests that they are small. Therefore, in 
what follows we will use the results for the form factors obtained in the static limit of QCD to 
compute the decay widths, but in the computation of the phase space we will use the physical 
meson masses. 

3 Result for 7 3 / 2 (l) is read from the plot in ref. [5], while the result for 71/2(1) was presented in ref. |14j . 
4 Results we quote in eq. JSJ are obtained after converting the values from ref. [15] to our definitions of Isgur-Wise 
functions, namely, 73/2(1) = r(l)/y/3 and 71/2(1) = C(l)/2, where 7(1) and £(1) are denned in ref. [16] . 
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Class I (a) Class I (b) 

Figure 1: Diagrams contributing to the Class I non-leptonic decay B — »■ D it: (a) pion emission, (b) weak 
annihilation. 



3.5.1 Semileptonic decays in theory 

The branching ratio of the semileptonic i?-decay to a j p = (l/2) + state should be very small 
compared to the decay to a j p = (3/2) + meson. A suppression due to the IW functions 

^^-017 (9) 
|r 3/2 (l)| 2 ' ( j 

is further enhanced by the phase space suppression (c.f. ref. [1]), and the suppression becomes 
one order of magnitude. Note that the decay to Di is less reliable because at w = 1 its amplitude 
is zero. 

Using the results of the quark model calculation in the BT formalism with the GI potential 
model one has [3]: 

B{B° d -»• D*~£u) ~ 0.7 x 1(T 2 , 
B{B° d -> D-, 3/2) £u) ~ 0.45 x 10" 2 , 

B(B° d -+ D'- {m lv) ~ 0.7 xHT 3 , 

B{B° d -»• D*-£u) ~ 0.6 x 10~ 3 . (10) 



Finite width effects are not negligible in the case of broad states, but they would reduce the 
predictions (by about 20%), thus further aggravating the problem we are addressing, i.e. the 
problem that predictions seem to be too small with respect to experiment. 

3.5.2 Non leptonic B — > D**tt + decays in theory 

Semileptonic decays would in principle provide the cleanest test of the theoretical predictions, 
but the undetected neutrino prevents us from doing a very good analysis. The above predictions 
can fortunately be tested by considering the non-leptonic decays, if an extra assumption is made, 
namely factorization. 
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u u 
Class III (c) 

Figure 2: Diagrams contributing to the Class III non-leptonic decay B — > D tt: (a) pion emission, (b) weak 
annihilation, (c) D -meson emission. 

As it is well known, there are three classes of non-leptonic decays. B^ — > D ' 7r , for 
example, belongs to Class I and is described by the sum of two diagrams: the pion emission 
through W (which is colour favored, c.f. fig. [Tk), and the annihilation through the W exchange 
(shown in fig. Q~b). The annihilation is expected to be small and the pion emission can be easily 
evaluated in the factorization approximation as a product of the annihilation constant f n and the 
B — > D form factor. As before, we use the form factors computed in the heavy quark limit, 
whereas in the phase space computation we use the physical meson masses. Using the values 
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given in eq. © one has [3] : El 



B(B° d -> D*-tt+) ~ 1.1 xHT 3 , 
£(B° ^ Z? M3/2) 7r+) ~ 1.3 xlO" 3 , 

B(B°^ j d;- (1/2) tt + ) ~ 1.1 xlO" 4 , 

B(B° d ^D*-ir + ) ~ 1.3 xlO' 4 , (11) 

where we include the ^-dependence of t 1 / 2 ,3/2( u; )) away from w = 1, which reduces the rate by 
around a factor of 2. The qualitative picture one gets from this exercise is that, similarly to the 
case of semileptonic decays, the decay rates to j = 1/2 states should be an order of magnitude 
smaller with respect to those with j = 3/2 in the final state. 

If one considers a class III decay, such as B + — > D ' 7r , then a priori three diagrams 
contribute: (i) pion emission through a color suppressed VK-exchange (see fig. [2^), (ii) annihilation 
of B through W, shown in fig. [2b, which is negligible because of the factor oc V u b, (iii) emission 
of the D meson through W-exchange (see fig. [2b). Although color suppressed, the last diagram 
cannot be neglected for the decay to j = 1/2 because its size is similar to the pion emission. This 
is a consequence of the smallness of ^^{w) [3]. On the other hand it vanishes for D 2 , where the 
factor /£>*, that appears in the factorized expression of the amplitude, vanishes because the 2 + 
state does not couple to the weak current. Notice that fr> 1 = f 3/2 is also expected to be small, 
based on the heavy quark symmetry. 

Since there is only one sizable contribution, class I decays should be preferred to test the 
theoretical estimate of the B — >■ D ' form factors. Class III nevertheless offers an additional 
qualitative test, because the additional diagram leads to a large difference with class I. In principle, 
the cleanest way to assess the magnitude of Ti/ 2 (w) and Ti/ 2 (w) would be through the study of 
semileptonic B — > D** decays, because the experimental extraction from the class I nonleptonic 
decays could be spoiled by the presence of the resonant nir pair in the final state. 



4 Experimental situation 



In contrast to the consistency of theoretical approaches, we find a rather different situation 
on the experimental side, especially in semileptonic decays where blatant inconsistencies are 
found between two sets of measurements by BaBar and Belle. We begin by explaining why 
semileptonic decays are in principle more difficult to analyze than non-leptonic three-body decays 
like B —7- D 7T7t, which may seem paradoxical since the former have a much larger rate. 

Dalitz plot analyses of B — > D -kit decay channels at B-factories have provided informations 
on the production rate and on the resonance parameters of broad D** resonances. Events are 
selected if the energy of the candidate is compatible with the beam energy and if the mass of 
the system formed by its decay products is compatible with the nominal B-meson mass. For an 
integrated luminosity of 500 fb _1 and an assumed decay branching fraction of 10 -3 there are 
typically 4000 and 9000 reconstructed signal events for the D** ->■ D*+n~, D* + ->■ D°tt + , D° -> 
K~ir + , K~ir + ir + Tr~ and D** — > D + ir~ , D + — > K~ir + ir + decay chains, respectively. 

Because of the missing neutrino, i?-meson semileptonic decays are more difficult to analyze. 
It is necessary to fully reconstruct the other .B-meson (Bt ag .) and a cut on the missing mass 

5 The expressions for the amplitudes involve the coefficient ai [T7], for which we take <zi ~ 1. 
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squared is used to select events with only a missing neutrino (Belle) or keeping also events in 
which the soft pion from the cascade D* — > Dir escapes detection, in addition to the neutrino 
(BaBar). These analyses have an efficiency which is typically two orders of magnitude lower than 
for the exclusive B — > 3-body decays considered previously. In practice, because semileptonic 
branching fractions into individual D** states are an order of magnitude higher than in exclusive 
non-leptonic final states, there is typically only an order of magnitude difference between the 
statistics of signal events analyzed in non leptonic and semileptonic S-meson decays. 

4.1 A long-standing confusion in semileptonic decays 

Unless explicitly stated, the numbers presented in this section are obtained by using the values 
given by the HFAG collaboration p~8|, and we average measurements from neutral and charged 
£>-meson using isospin symmetry. Obtained values are quoted for the B d meson. 

The inclusive semileptonic decay branching fraction of B d and of B + decay is far from being 
saturated by the sum of D£ + U£ and D £ + vi decay channels. More specifically |19| : 

B{B^XJ + u e ) = (10.09 ±0.22)%, (12) 
B{B%^Dl + v e ) = (2.12 ±0.06)%, 
B{B° d ^D*l + u t ) = (5.11 ±0.10)%. 

In other words, the semileptonic branching fraction to the charm states which are not simply a 
D or a D* is thus equal to: 

B(B° d -> non-D { * ] £ + v e ) = (2.86 ± 0.25)%. (13) 
Decays to the narrow D** states have been measured with good accuracy: 

B(B° d -> D* 2 £ + u e ) = (0.29 ± 0.03)% , 

B(ti%->Di£+v t ) = (0.58 ±0.05)%, (14) 

giving: 

B(B° d -> DZrroJ + ve) = (0-87 ± 0.06)% . (15) 
The above values include the branching fraction of D** into the observed final state. H Another 
piece of information comes from the measurements of the exclusive B — > D ttVv^ decays: 

B(B%-*ThPv t ) = (0.60 ±0.06)%, 
B{B%^TTirt+vt) = (0.83 ±0.06)%, (17) 

giving: 

B(B° d -> D^iri+Uf,) = (1.43 ± 0.08)%. (18) 

This value can be compared with the expected D and Dtttt£ + i'£ branching fractions from 

the decays of narrow D states given in Table [TJ 

6 Few branching fractions of D** decays into exclusive final states are not well determined and we use the 
following values: 

B(D* 2 ° -> D + tv~) = 0.41 ± 0.02 , B{D? -> D* + tv~) = 0.26 ± 0.02 , B{D* ° -> D + tt") = 2/3 , 
B(_D? -> D* + tt~) = 0.45 ±0.02, Z3(_D?' -»■ D* + tt~) = 2/3 . (16) 

To make these evaluations we have assumed in addition that: (i) decays exclusively into D-k or D*tt. Channels 
with two charged pions have been studied and no signal was observed; (ii) Do decays exclusively into Dir; (iii) Di 
decays into Dmv and D*tt with a ratio B(D" — > D a ir + iT~)/B(Di — > D* + tv~) = 0.32 ± 0.03, and we assume that 
the decay proceeds through the chain Di — > Dq-k. 
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Decay channel 


Di 




total 


B{B u d -> Diri+ve) 




(0.18 ±0.02)% 


(0.18 ±0.02)% 


B(B° d -> S'tt^+i//) 


(0.39 ± 0.04)% 


(0.11 ±0.01)% 


(0.50 ± 0.04)% 




(0.19 ±0.02)% 


0.00 


(0.19 ±0.02)% 



Table 1: Branching fractions for B® — > D ■n{n)i + vi decay channels where the hadrons cascade from a 
narrow D meson. 



From these measurements one can draw several conclusions: 

• narrow D\ and D\ states, with no additional pion, account for about 1 /3 of the non-D^ £ + V£ 
final states; 

u\ u\ 

• D tt£ + V£ final states account for about 1/2 of the non-L> l + vg final states. As a result, 

final states with two or more pions should account for the other half; 

• the broad state component of the Dir system corresponds to a branching fraction equal to 
B(B° d -> [Dir] broad £+u e ) = (0.42 ± 0.06)%; 



• the broad state component of the D ir system corresponds to a branching fraction equal to 
B(B° d -> [D\] broad £+u e ) = (0.33 ± 0.07)%; 

u\ 

For theorists, it remains to interpret the origin of the broad D nir components, with n > 1 
which correspond to 2/3 of these hadronic final states in B semileptonic decays. 

There is, at present, an apparent contradiction between the measured values for the D Q . 

B(B° d ^D* £ + ^) = (0.35 ±0.07)%, (19) 

and the corresponding theoretical expectations. According to theory, the production of these 
broad resonances should be much lower than the one of narrow states and this is apparently not 
verified (see below for details). For the broad D 1 state, the situation is different because the two 
experiments disagree. Belle does not see any broad D v component, while BaBar gives : 

B{B Q d -> D'^+ve) = (0.26 ± 0.04 ± 0.04)% , 

HFAG gives ("Updates of Semileptonic Results for End Of 2009") 

B{B° d -> D' 1 £ + u e ) = (0.13 ± 0.06)% , 

but it must be understood that the two measurements are incompatible (BaBar and Belle results 
differ by 3.2 a.) The PDG group discards Belle without explanation. 
Meanwhile several comments are in order: 

• experimenters cannot claim that they have really measured the production of the broad D*q 
and Di resonances. There could be additional contributions from broad Dn and D*7r final 
states in the registered spectra; BaBar states explicitly that they have not subtracted any 
non resonant background, for lack of a satisfactory fit for it; 

• the branching fraction attributed to the Dq is compatible with the broad component rate 
obtained by analyzing the Dir final state; 
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• for the Di production, the quoted value of BaBar is compatible with the broad component 
rate obtained by analyzing the D ir final state. 

All this is compatible with the idea that the real difficulty causing the disagreement within 
experiment, and perhaps with theory is the one of analyzing events in terms of broad resonances, 
as we discuss in subsection 15.41 

4.1.1 Summary 

In Table [2] are summarized the present measurements of i?-meson semileptonic decays into a 
charm hadronic system. Values are quoted for the B d ; corresponding results for the B + can be 
obtained by multiplying these values by the lifetime ratio t(B + )/t(B®) = 1.079 ± 0.007. 



Decay channel 


branching fraction (%) 


B" d -> X c £+u e 


10.09 ±0.22 


B° d -»• D-i+ui 


2.12 ±0.06 


B° d -> D*-£+v e 


5.11 ±0.10 


B° d -> D^t+vi 


0.58 ± 0.05 


B° d -> D^-e+vt 


0.29 ± 0.03 


B° d -> [Dir] broad l+v t 


0.42 ± 0.06 


B° d -+[D*n] broad e + v t 


0.33 ± 0.07 


B° d -> [D7nr] narrow l + vi 


0.19 ±0.02 


B° d -+ x~r^ 

B° d K°£ + V£ 


1.24 ±0.26 


0.06 ±0.01 [2D] 



Table 2: Measured semileptonic B d branching fractions. The [£)7r7r] narrow hadronic final state corresponds 

to the decay of the D±. The ^^broad" 18 hadronic final state contains a D or D meson with at least two 
pions or aijora rf meson. 

From these measurements there are at least two questions which remain to be clarified: 
u\ 

• the origin of [D Abroad states. What is the fraction of these states which can come from 
the D*q and D l mesons? A possible answer to this question is the subject of the present 
paper. 

• the contribution of broad final states with several pions or with a r/ or a r( . Because of the 
large mass of the rj- 1 ^ mesons it is not expected that corresponding final states have a large 
contribution. 

4.2 B ->■ D**tt + decays 

In this subsection we provide a summary of present measurements at BaBar and Belle of the 
decays B — )■ D tt + . 

BaBar and Belle collaborations have measured several B — > D ir + decay channels using 
Dalitz analyses. Averaged values of B — > D n branching fractions measured by BaBar [21|, [22] 
and Belle (231 EU [25] are given in Tabled 

A few remarks can be made: 
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Decay channel 




B+ 


W 2 n+ 


(4.9 ± 0.7) x 10~ 4 


(8.2 ± 1.1) x 10" 4 


D l7 r+ 


(8.2+^) x 10~ 4 


(15.1 ±3.4) x 10~ 4 




< 1 x 10™ 4 


(7.5 ± 1.7) x 10" 4 


w n+ 


(1.0 ±0.5) x 10" 4 


(9.6 ±2.7) x 10" 4 



Table 3: Measured branching fractions for B — > D n + decay channels. 



• branching fractions are higher for the B + than for the B®, where both are measured, Lj 

• considering the D 2 production, which is the most accurate, it is also not too far from 
equality, as would be expected according to factorization, since there is no diagram with 
D 2 emission. On the contrary, it is expected that for the + the two rates should be very 
different, as it is indeed found (see below). 

• D\ production seems to be higher than D 2 , in a certain contradiction with heavy quark 
symmetry. This is understandable by a simple l/m c effect, as in semileptonic decays. 

• the production of Dq states is not well measured. In B + decays it seems to be similar to 
the D 2 but in B® decays it seems to be much smaller. In fact, measurements of B® decays 
from Belle and BaBar (preliminary) are not in good agreement even if it is difficult to draw 
a clear conclusion because of the attached uncertainties: 

B{B° d -»• D* ~tt + ) x B(D*~ -»• D°tt~) = (20) 
(0.60 ± 0.13 ± 0.15 ± 0.22) x 10~ 4 Belle , 
(2.18 ± 0.23 ± 0.33 ± 1.15 ± 0.03) x 10~ 4 BaBar . 

BaBar reports a larger systematic uncertainty, coming from the modelling of the fitted 
distribution, than Belle. Anyway, the decay of neutral B® is in both experiments clearly 
smaller than the charged one, and this can be understood theoretically because in the 
charged case, and, contrarily to D 2 *' , there is a diagram with emission of Dq*' which can 
overwhelm the pion emission diagram, which is small because of the smallness of Ti/ 2 (1). 



5 Comparison between theory and experiment 

Results of the preceding sections are summarized in Table 01 Let us then recapitulate the con- 
clusion one can draw by taking the experimental data as they are presented. 

5.1 Ratio of B° d ->■ D*-£ + u e and B° -> D*-n + 

Assuming the validity of the QCD factorization and by describing the B — > D* transition matrix 
elements by a slowly varying Ti^^/^iw), one can easily see that B d — > Dq~tt + and B® — > Dq~£ + V£ 
decays are governed by T\i 2 alone. |£| Using the values given in Table [H the ratio of semileptonic 

7 In these comparisons between branching fractions for charged and neutral B-mesons we are interested in 
differences which appear in addition to the 7 % expected from the lifetime difference. 

8 The general idea of the relation between semileptonic and non leptonic decays is due to M. Neubert [26] . 
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$3theovy 




^expt. 1 ^theovy 


B°, -> TT"e + v e 








A 


0.7 x 10~ 2 


(0.29 ±0.03) x 10" 2 


~ 0.5 


Di 


0.45 x 10~ 2 


(0.58 ±0.05) x 10~ 2 


- 1. 


d\ 


0.7 x 10~ 3 


[0..3.2] x 10~ 3 


[0.,5.] 




0.6 x 10~ 3 


(3.5 ±0.7) x 10" 3 


6. ±1. 










Tf 2 


1.1 x 10~ 3 


(0.49 ±0.07) x 10~ 3 


- 0.5 




1.3 x 10~ 3 


(8.2±f;f) x 10- 4 


[0.5,1.] 




1.1 x 10~ 4 


< 10- 4 (90% C.L.) 


no result 


D* 


1.3 x 10~ 4 


[0.3,3.4] x 10^4 


[0.2,2.6] 



Table 4: In this Table are collected the values expected and measured for D production in semileptonic 
and non leptonic B® meson decays. These values have been given already in previous sections. The 
theoretical expectation is taken to be the one of the quark model, subsection 13.51 A range of values is 
given within brackets when there is not a good compatibility between BaBar and Belle measurements. In 
this case we take the minimum value minus one sigma and the maximum value plus one sigma to define this 
range. In general there is agreement between measured and expected branching fractions for narrow states. 
For broad states results are in contradiction with expectations (mainly the D production in semileptonic 
decays) or rather uncertain. 



to non-leptonic decays with D in the final state must be ~ 5. Experimentally, instead, such a 
ratio spans a large interval between 8 and 140. In contrast to that situation, decays to the narrow 
D 2 state lead to a ratio that is theoretically expected to be equal to 6, which is confirmed by the 
experimentally established value 6±1. For decays to D\ state, uncertainties are larger and based 
on a single unpublished result from Belle but the expected theoretical value for the ratio, which 
is equal to 3.5, agrees roughly with experiment (7 ± 2). 

5.2 Contradiction between the phenomenological predictions and the semilep- 
tonic experimental data 

Now, we can go further still and state that the semileptonic experimental data contradicts the 
HQET estimate for the decay to a j = 1/2 state, with a huge discrepancy which is one order of 
magnitude in rate. 

To arrive to such a conclusion, one first has to take into account the disagreement among 
experiments in B(B — > D l i + V() states. While the result reported by Belle seems to be compatible 
with expectation of a very small rate, the result of BaBar is much larger and disagrees with both 
Belle and the expected value. Both experiments, instead, agree on the value for B(B — > D £ + vi) 
which is far too large when compared with expectations. While the results by BaBar are far too 
large when compared to the expectations, they are still consistent with the heavy quark symmetry 
expectations, i.e. the two rates are nearly equal. The results by Belle instead indicate a complete 
breakdown of the heavy quark symmetry. On the whole, it is fair to say that both experiments 
disagree with theory for both j = 1/2 states. 

On the other hand, there is a qualitative agreement in both types of transitions to j = 3/2 
states. There is an excess of theory, by a factor two, for B{B — > Z? 2 ^ + ^), but there is also an 
overall success for the sum £(£^[z4DiK+^)~l%. 
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5.3 Better situation for non leptonic decays, yet not conclusive 

The situation with non-leptonic decay to a j = 1/2 state is much better not only in experiment, 
but also concerning the comparison between theory and experiment. For the Class I decay, 
B^ — > Dq~7t + , the prediction pip coincides with the Belle measurement, and is compatible 
with BaBar within the quoted uncertainties. Notice the important point that in Class I decays 
factorization is expected to hold to a good approximation both on theoretical grounds and also, 
taking into account a large number of decays with such topology, on empirical grounds. 

The discrepancy between Belle and BaBar occurs in B® — > Dq~tt + , which could be attributed 
to the difficulty of extracting a broad resonance, with possible large non-resonant structure, and 
with the additional difficulty of a irir crossed channel interference (see below). 

A fact that seems to attest the soundness of the theoretical statements about the smallness of 
the production of the j = 1/2 states is the large difference between neutral and charged -B-decay 
into the broad D -state: the charged decay rate is much larger than the neutral one, by about one 
order of magnitude as given in Table [3] This is easily understood because an additional diagram 
is present in the charged case, which is the D -emission (class III) - El Although color suppressed, 
this diagram gives a contribution much larger than the one with the pion emission, if Ti/2(1) is 
small [3]. In that case the D** emission amplitude dominates the charged rate, and dominates 
over the neutral decay amplitude. A similar effect is observed in the case of the broad 1 + final 
meson. Although a branching ratio has not been published, the bound on the neutral -B-decay 
in ref. [27] clearly indicates that the charged decay is much larger than the neutral one. 

The discrepancy of around a factor of two between charged and neutral B decay to D 2 , could 
be interpreted as an estimate of the correction to the factorization approximation in which the 
two decays are expected to have nearly equal rates. Such a discrepancy is similar to what is found 
in common tests of factorization [3j CE7] EHJ [29l SO] (see also references therein) . 

5.4 Discussion of the main discrepancy and possible explanations 

If we believe the results of theory, which are rather consistent, and if we take the experimental 
results for broad states in semileptonic decays, then one or both states have much too large 
rates as compared to theory. One experiment also suggests a complete breaking of heavy quark 
symmetry. In non-leptonic decays there is a better agreement between theory and experiment 
but present uncertainties in B^ — > D it decays are too large to derive firm conclusions. 

Of course, one could evoke weaknesses in the assumptions which allow to derive phenomeno- 
logical predictions. In particular one can argue that the l/m c effects could be large. However, 
large l/m c effects cannot explain the contrast between a relative success in non-leptonic decays 
where they should be present too. One could also complain about he validity of the factorization 
approximation, but that is unlikely to be the case as factorization in the Class I decays passed 
many experimental tests and no large deviations have been found so far. Finally, let us stress the 
satisfactory qualitative agreement in the case of decay to a j = 3/2 state, both semileptonic and 
non-leptonic ones. 

The problem of broad resonances . A possible reason for the qualitative agreement be- 
tween theory and experiment in the f?-decays to a j = 3/2 state can be explained by the fact 
that the j = 3/2 states are narrow. Distinguishing very broad resonances from continuum is 
extremely difficult enterprise, both on theoretical and experimental sides. 

There is no unambiguous way of writing the broad resonance line shape, all the more for 

9 Such an explanation was first offered by Belle [27] • 
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S-wave scattering where very strong couplings can be present, and therefore the very notion of 
separating a resonance and the non-resonant continuum is theoretically ambiguous. Furthermore, 
the qq states could be competing with non qq states in S waves and additional resonances could be 
generated by the scattering. Finally one can also encounter problems with contributions arising 
from the tails of the ground state (denoted as D v , B* in tab. II of ref. [25]) or of radial excitations 
in D { *\. 

Ideally, one should be able to compare the whole amplitude with experiment, and not just 
the resonance under study, but that is obviously not possible in practice. All this underlines the 
advantage of working with narrow resonances. 

It must be repeated, however, that if broadness was a sole cause for a large discrepancy 
discussed above, then one would be short of explanation regarding the non-leptonic decays for 
which the disagreement is not large. Keep in mind, however, that potentially large uncertainties 
due to the arbitrariness of the non-resonant continuum should enter the game also in the non- 
leptonic case. Last but not least, for neutral — > D^~tt + — > D°7r~ir + , which is the relevant 
one for our purpose, one can have interference with the crossed channel Trn which resonates into 
p,fo,etc... (Bj — > D p°, ...) All these contributions cannot be separated out without heavily 
relying on specific models and the resulting uncertainty may lead to inconclusive comparison 
between theory and experiment. 

Blaming broadness of states for the difficulties in measuring the rates of j = 1/2 is strongly 
supported by the following argument. In B® — > D 7r~7r + , Belle and BaBar find exactly the 
same total rate, and the same rate for all the decays to relatively narrow resonances, i.e. not 
only B — > D 2 tt, but also B — > Dp, B — > D/2(1235). On the other hand, large discrepancies 
appear in the central values of the decays to broad resonances, not only in B — > D ir, but also in 
B Df (600) (S-wave). 



Part II 

Proposal for the complementary study of the narrow 
strange counterparts 



6 Motivation 

Our proposal starts from the above observation that analysis of broad resonances has always been 
a difficult task. The fact that no special problem arises for the narrow j = 3/2 states suggests 
that the broadness of j = 1/2 states in the non-strange case could be the origin of the difficulties. 
At least, it could help much if one could deal with states analogous to the controversial D** (i.e. 
Dq and D[), but narrow. Even if not leading to an immediate solution, it would substantially 
help clarifying the comparison between theory and experiment. 

Furthermore, a study of B® — > D*^~7r + would be an important test of the consistency between 
theory and experiment as far as t 3 / 2 (1) is concerned. 

10 Note that PDG use the notation /o(600) for the lowest scalar J PC — ++ state [31], that is often referred to 
as cr(600), or e(600). 
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6.1 The two narrow j — 1/2 D s j states 

It is very fortunate that the strange j = 1/2 D**-states, -D* (2317) and D s i(2460), are very 
narrow, because their masses are below their respective D^K thresholds. The broad non-strange 
states are heavier than the D^tt threshold. While the SU(3) symmetry breaking is large in the 
phase space, it can still be expected to work well for the electroweak amplitudes and strong 
couplings, as has been observed most often. 

The narrowness of the states offers an exceptional possibility to test the theoretical predictions 
in a much better experimental situation. It eliminates at the same time the problem of the non- 
resonant background, and interference with competing crossed channels, since both should be 
relatively negligible near the peak. 

The effect of SU (3) breaking is expected to be small for the lowest lying states with given 
quantum numbers. We therefore expect Ti/2(1) to be rather close to the non-strange case. Note 
that in the lattice QCD study of ref. [TT], no significant dependence of ^^(l) on the light quark 
mass has been observed. 1^1 

We should emphasize once again a great advantage of the non-leptonic over semileptonic B®- 
decays in that they do not have the neutrino identification problem, but have the two body final 
state with well known masses. Theoretically, B® — > D s j it is the most interesting decay because 
it is described by the pion emission diagram only (B s annihilation being neglected as usual). In 
the factorization approximation, it directly yields Ti/2(1)- 

Warning concerning a possible misinterpretation of D s (2317, 2460): A potential caveat con- 
cerning the D* (2317) and -D s i(2460) is that they might not be the ll qq" states. A controversy 
resides in the fact that the measured masses of these states are lower than predicted. However, 
the level ordering of the ll qq" states, 0~, 1~ , + , 1 + , 1 + , 2 + , is consistent with what is observed 
with the D s j mesons so far. Moreover, the study of their transition properties does not favor an 
exotic assignment either. We must underline that a measurement of the decays proposed here 
will also provide an extra check of the qq structure of D s (2317, 2460). 

7 Decay branching fractions of D* (2317) + and D sl (2460) + states 

Of course, to measure the B® — > D~^j7r~ rates, the knowledge of the D s j branching ratios is 
necessary. 

In ref. [31] only absolute values for the -D^(2460) branching fractions are quoted. This is 
because, at present, there is only a single measurement [M] of the D^(2460) production in 

B — > 2)^(2460)1) ^ decays, independently of the decay channel for the D^(2460). Production 

of D s j states was studied by considering the missing mass distribution in B — > D X decays and 
signals were observed only for X = Df, D* + and 1)^(2460). As a result, there is no absolute 
decay branching fraction measurement for the Z?* (2317). 

7.1 D s * (2317) + decay channels 

Experimental results collected in [31] are reminded in Table [5j 

11 A proposal to study the B s — > D s j transition has been made in ref. [32] in order to test whether or not the 
D 3 j states are indeed the u qq" structures. They use the QCD sum rule calculations in HQET and find huge SU(3) 
breaking effect (~ 100%) in the form factor [compare eq. (34) in ref. [15] with eq.(32) in ref. [32]], which contradicts 
the lattice QCD findings of ref. [11] , 
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Decay channel 


90% C.L. limit 


comment 


Dtl 


< 0.05 


forbidden 




< 0.059 


allowed 


L>+77 


< 0.18 


allowed 




< 0.11 


forbidden 


D+7T+7T- 


< 0.004 


forbidden 


D+ir°TT° 


< 0.25 


forbidden 



Table 5: 90% C.L. limits on branching fractions for different decay channels measured relatively to the 
Dfir channel. In the last column are indicated the allowed and forbiddden decay channels from angular 
momentum and parity conservation. 



The electromagnetic Df'yy is expected to be negligible as two photons have to be radiated. 
It thus remains only two possible decay channels for the -D* (2317) + . In Table [6] are indicated 
some model expectations on these decay channels [351 136] . 



Decay channel 


model 1 1351 


model 2 [36] 


90% C.L. limit 


D* s + 1 


92.5% 
7.5% 


84% 
16% 


< 0.059 



Table 6: Some model expectations for _D* (2317) + branching fractions compared with the experimental 
result. 

The present limit on the D* +, y decay channel is more stringent than the estimates; in the 
following we will use: 

£(£>* (2317)+ -> D+ir°) = (97 ± 3)%. (21) 
7.2 D s i(2460) + decay channels 

Experimental results collected in [51] are reminded in Table 



Decay channel 


value or limit 


comment 


Dtl 
D s * (2317)+ 7 

D+TT° 
D+7r 7T 

D+jj 


(48 ±11)% 
(18 ± 4)% 
(4.3 ± 1.3)% 

< 0.08 (90% C.L.) 

(3.7±|;2)% 

< 0.042 (95% C.L.) 

< 0.68 (95% C.L.) 

< 0.33 (95% C.L.) 


allowed 
allowed 
allowed 
allowed 
allowed 
forbidden 
allowed 
allowed 



Table 7: Measured branching fractions or upper limits for different D s i(2460) + decay channels. In the 
last column are indicated the allowed and forbiddden decay channels from angular momentum and parity 
conservation. 

Many decay channels are possible and individual decay branching fractions are not accurately 
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measured. The situation is thus experimentally less favourable than for the Z?* (2317) + resonance 
to measure the production rate of this state. 

8 Expected rates at LHCb 

8.1 Analysis method 

We would like to have a measurement of the decay chain B® — > D*^7r + , — > D~ir°, D~ — > 
K + K~ir~ in which the 7r° meson cascading from the D*^ is not detected. 

It is proposed to measure the missing 7r° 4-momentum using the measurement of the B® 
direction and two mass constraints (m^o and m B o). The B® direction is determined from the 
reconstructed positions of the pp interaction and the B® decay vertices. Measured uncertainties 
on these quantities can be included in a fit with the two mass constraints. 

There could be 2 solutions for the signal and a study based on simulated events may help to 
choose one of these possibilities. The amount of background candidates can be decreased using 
the fit x 2 probability. 

For signal events, as the D*q has a very small intrinsic width, one expects to observe a peak 
in the D~tt° mass distribution having a width which depends mainly on the accuracy of tracking 
capabilities. 

8.2 Expected rates 

The proposed analysis is based on the same charged particles final state which was already 
measured in LHCb for the channel: B® — > D~ir + , D~ — > K + K~ir~ . Few selection criteria have 
to be removed to allow for the missing 7r° meson and in particular the condition on the similarity 
between the directions defined by the two vertices and by the K + K~tt + it~ momentum. 

Analyzing 336 pb~ l integrated luminosity, LHCb has measured [37] about 6000 B® — > D~tt + 
decays. The number of B® — > D*q1t + reconstructed events can be estimated by comparing the 
corresponding branching fractions for the two decay channels. 

From SU(3) symmetry and factorization, we can simply identify the branching fraction of 
B s — > D s jtt with the one of the neutral B into charged D and tt. Indeed, the phase space is 
also very close to the one in the non strange case. In view of the other uncertainties, we can 
safely disregard any SU(3) effect. This means from the measured case, the one of J p = + , 

B(B° S -> D*^(2317)tt + ) = (1.0 ± 0.5) x 10" 4 , (22) 

where we average the results of Belle and BaBar for the non-strange decays (BaBar is presently 
not published). This value agrees with the theoretical expectation using the heavy quark limit 
[10" 4 , c.f. eq. ([TT]) ]. However, using the experimental value for the non-strange decays together 
with the SU (3) light flavor symmetry is likely to be better than the result derived in the heavy 
quark limit and assuming exact factorization. 

To assess the soundness of the SU(3) assumption, let us consider the decays to D,D S . The 
LHCb collaboration has measured: 

B{B° S D~ir + ) = (2.95 ± 0.28) x 10" 3 . (23) 

In this expression we have added in quadrature the different uncertainties quoted in the pub- 
lication. The value agrees well, as expected, with the corresponding measurement for the B^ 
meson. 

B{B° d D~tt + ) = (2.68 ± 0.13) x 10" 3 . (24) 



18 



Analyzing an integrated luminosity of 1/6 , the LHCb collaboration can thus expect to 
reconstruct: 

N(B° a -»• L>^(2317)vr + ) = 1800 xjx(l± 1/2) x B(D* s0 -> L^vr ) x e w o, (25) 

with the Z?~ meson reconstructed in the K + K~n~ decay channel. The quantity e^o corresponds 
to the efficiency of the additional cuts which have to be applied to select the events. 

A very few hundred of events are expected and the signal visibility will thus depend mainly 
on the mass resolution for the D~tt° system and on the combinatorial background level. 

9 Conclusion 

9.1 Feasibility of the proposal 

We propose an experimental study of the B s — > D s jir decays that would provide us with an impor- 
tant verification of the observations made in the corresponding non-strange modes. Furthermore 
it would allow us to elucidate the problem of small value of 7i/ 2 (l). 

If a really unexpected value for B(B® — > D*q-k + ) is found, this could mean that 

- either we are mistaken in the theoretical evaluation of 71/2(1); which would be very surpris- 
ing in view of good consistency of several approaches, or the l/m c corrections are exceedingly 
large in the j = 1/2 case, 

- or the narrow D s j-stdXes situated below the D^K thresholds are not the "gq" states with 
j = 1/2 (ree ref. [38] for a review). 

Both these possibilities do not seem plausible to us. The remaining uncertainty on the theoretical 
side could be significantly reduced by the lattice study of the B® — > D s transition form factors 
at finite heavy quark masses. 

If the expected rate is confirmed, that would set beyond doubt the theoretical estimates of 
small values for 7i/ 2 (l) and it would confirm the assignment of the D s j states. A strong suspi- 
cion would be confirmed against the semileptonic measurements or identifications of resonances 
performed in the non-strange case. 

9.2 Remaining problems on the non-strange side 

Even if the answer of the proposed experiment is in agreement with theoretical expectations made 
by adopting the "qq v assignment to the D** states, it will still not give us the full explanation to 
the problems observed in the non-strange case. The problems encountered on the experimental 
side, especially in semileptonic non-strange decays, remain to be understood: the origin of the 
discrepancy between Belle and BaBar; why so large apparent rates for decay to + ? A theoretical 
explanation for the large number of events in the non-strange semileptonic decay is missing. 

The observed excess of events in D^tt (around 1%) and in Dnir, that in our opinion are 
not the lowest j = 1/2 or j = 3/2 states, needs an explanation. Such events should have their 
counterpart in non-leptonic decays. To test an excess in the Dn channel, a study of the decay 
B® — > D°7r~7r + at LHCb would be very welcome. 

12 To interpret these events, one could think of a possible contribution from radial excitations in the non-strange 
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